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ABSTRACT 

We present a self-consistent, relativistic model of rapidly rotating neutron stars de- 
scribing their exterior gravitational field. This is achieved by matching the new solution 
of Einstein's field equations found by Manko et al. ( 2000) and the numerical results for 
the interior of neutron stars with different equations of state calculated by Cook et al. 
(1994). This matching process gives constraints for the choice of the five parameters 
of the vacuum solution. Then we investigate some properties of the gravitational field 
of rapidly rotating neutron stars with these fixed parameters. 
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1 INTRODUCTION 

Accretion on to neutron stars is primarily determined by 
their exterior gravitational field for which an analytical de- 
scription is preferable to a numerical one. The solutions of 
Einstein's field equations suitable for Black Holes as the 
Schwarzschild- and the Kerr-metric do not hold for neutron 
stars because those have a solid surface and are therefore 
oblate when rapidly rotating. The second difficulty is the 
possible existence of an intrinsic magnetic field. A new solu- 
tion has been strongly needed including the two parameters 
of the Kerr-metric, mass and angular momentum, and three 
additional ones, the mass quadrupole moment due to the 
oblateness of the star, the electric charge (which is of no im- 
portance for astrophysical objects) and the magnetic dipole 
moment. After more than ten years of work Manko et al. 
( 2000 ) found a space-time meeting all criteria. 

These parameters are in a mathematical sense for- 
mally independent, but in an astrophysical sense, the mass 
quadrupole moment is a function of the mass and the angu- 
lar momentum. To find these dependencies we use numerical 
results found by Cook et al. 11994t . They solved the general 
relativistic structure equations for an axisymmetric rotating 
star for 14 different equations of state (EOS) and three se- 
quences and presented certain quantities for the stars in a 
tabulated form for five of them. By comparing those with 
the analytical values we get correlations for the parameters 
of the new solution. 

This was also done by Sibgatullin & Sunyaev (1998), 
but with another metric IManko et al. 1994t and only three 
EOS with low stiffness. These authors investigated their 
found models in further publications (the most recent is Sib- 
gatullin & Sunyaev (2000)). This metric, however, did not 
permit the rational function representation - the reason why 



Manko et al. continued their search for a new solution - and 
could not be written in a concise form which is more useful 
for analytical examinations. Furthermore the same proce- 
dure - but in opposite order - was done for this new metric 
extended to all five EOS and all three sequences. 

With fixed parameters this new solution is highly useful 
to improve the models of rapidly rotating neutron stars e.g. 
in accretion disc calculations for which the Schwarzschild- or 
Kerr-metric are still used at this time. These introduce large 
errors because of their significant deviations to real neutron 
stars due to symmetry simplifications which are now not 
necessary anymore. 



2 THE GENERAL EXTERIOR FIELD OF 
ROTATING NEUTRON STARS 

Starting with the most general form of a stationary axi- 
symmetric metric, found by Papapetrou (1953), 

ds 2 = / (dt -wd(ff - / _1 [e 27 (dp 2 + dz 2 ) + p 2 d<j>% (1) 

where p and z are quasi-cylindrical Weyl-Papapetrou- 
coordinates and /, w and 7 are functions independent of cj> 
and t, Ernst (©68 ) found that Einstein's equations reduce 
to the form 

(Re£ + |$| 2 )V 2 £ = (V£ + 2$*V$)- V£ 

(Re£ + |$| 2 )V 2 $ = (V£: + 2$*V$)- V$, (2) 

after introducing the two complex potentials £ and $ and 
their complex conjugates £* and $* as 

£ = if - \$\ 2 ) + if $ = A t + i<p', (3) 

(p and cp' are defined by 

/9 _1 / 2 Vw — 2h x Im($* V4>) = hxVip 
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p^fiVA^-wVAt) 



h X V<p' 



(4) 



with the t- and ^-components A t and A^, of the electromag- 
netic four-potential and the unit-vector in azimuthal direc- 
tion h. If the Ernst-potentials are known, the functions /, 
w and 7 can be directly calculated from 

/ = Re£ + $$* 

w, P = -AT 2 k z + 2/m($*$, z )] 
W, z = p/" 2 [v3,p + 2/m($*$, p )] 

7,p = ^r 2 [(f, P + 2$*$, p )(r p + 2$$: P ) 

-(£,« + 2 ***,,)(£,*,+ 2 
7* = Jp/~ 2 i?e[(£, p + 2<r$, p )(r z + 2$$;,)] 



-2p/- 1 J Re(c&; p $„ 



(•») 



The two Ernst-equations are solved by simple integral equa- 
tions ( Sibgatull in 199l| > which introduce the axis data of the 
Ernst-potentials e(z) = £ (p = 0, z) and f(z) = <!>(p = 0, z): 



1 / e(£)nid(T 

TV 



$ = 



(6) 



where pi has to satisfy 



(gW+e(C) + 2/(77)/(C))^ 

(C - n) 

1 

/" pi(a) dt 

J VT^ 



= 



da 



(7) 



with ?7 = z + ipr, £ — z + ipa, e(£) = e(|), r, cr G [— 1, 1]. 
In the first integral the Cauchy principal value has to be 
calculated. These new functions can be expressed as a mul- 
tipole expansion of the mass-, angular momentum-, electric 
and magnetic field-distribution 



£ = 



$ = 



(8) 



(n+1) 



C(p = o) 
9(P = 0) = 



(9) 



The physical multipole moments are the real (mass multi- 
poles) and imaginary (angular momentum multipoles) part 
of m n and the real (electric multipoles) and imaginary (mag- 
netic multipoles) part of q n . 

To derive a new exact analytical solution of Einstein's 
equations one starts with an ansatz of e(z) and f(z), checks 
whether the solution has the wanted parametrisation of 
the multipole moments and performs the cumbersome but 
straight-forward calculation of the metric functions /, w and 
7- 

Different authors (e.g. Aguirregabaria et al. (1993), 
Manko et al. 1199 4-)'! made an effort to find a new solution in- 
volving five parameters representing mass, angular momen- 
tum, charge, magnetic dipole moment and mass quadrupole 



moment of the star which can be expressed with relatively 
simple rational functions. After more than ten years and sev- 
eral ansatzes such a solutions was found (IManko et al. 20001 
starting with 



e(z) 

/(*) 
with 
8 = 

d = 



(z — m — ia) (z + ib) + d — 5 — ab 
(z + m — ia) (z + ib) + d — S — ab 

Qz + ip 

(z + m — ia) (z + ib) + d — S — ab 



(10) 



/' 



■2 1,2 



-(a- b) 2 - Q 2 
i 2 -(a-bf -Q 2 ] 



1 r 



(11) 



After checking the multipole moments one finds that m is 
the mass, a the specific angular momentum, Q the charge 
and p the magnetic dipole moment (if Q = 0). The param- 
eter b is no physical quantity, but is related to the mass 
quadrupole moment. For comparison the functions e(z) and 
f(z) of the Schwarzschild- and the Kerr-solution are 

-"' f(z) = (12) 



e(z) = 
<z) = 



z — m 



z + m z + m 

z — m + i a 



f(z) = 0. 



(13) 



z + m + ia 

To be able to write the metric function in a rational form it 
is convenient to introduce general spheroidal coordinates by 



V 



(14) 



with r± = p 2 + (z± k) 2 and k = ^/d + 5. The metric then 
becomes 



ds — f (dt — w 
dy 2 



2 r-1 



(x 2 



dx 
'\x 2 - 1 



+- 



i - y 

and the functions are 



) + (X 2 - 1)(1 - y W 



(15) 



J D 



with 

E = 



o 2 T 



E 



16k, s (x 2 



-(i-y 2 )f 

E 



(16) 



{4[k 2 (x 2 -l) + 5(l-y 2 )] 2 + (a-b) 
[(a - b)(d - S) - m 2 b + Qp](l - y 2 ) 2 } 2 
— 16k 2 (x 2 — 1)(1 — y 2 ){(a — 6)[«: 2 (x 2 — y 2 ) 
+25y 2 ] + (m 2 b-Qp)y 2 } 2 



D = {4(/t 2 



Sy 2 ) 2 + 2«;ma;[2K 2 (a; 2 — 1) 



+ (25 + ab- b 2 )(l - y 2 )] + (a- b)[(a - b)(d - 
-m 2 b + Qp](y 4 - 1) - 4d 2 } 2 + 4y 2 {2«: 2 (x 2 - 



■5) 
1) 



F = 



[Kx(a — b) - rnb] - 2mb8(l - y 2 ) + [(a - b) 

(d - 8) - m 2 b + Qp\(2nx + m)(l - y 2 )} 2 

8k(x 2 - l){(a - b)[K 2 (x 2 - y 2 ) + 28y 2 \ 

+y (m b — Qp)}{Kmx[(2nx + rn) — 2y (25 + ab 

-b 2 ) -a 2 + b 2 - Q 2 ] - 2k 2 Q 2 x 2 - 2y 2 (A5d 

-m 2 b 2 )} + {A[k 2 (x 2 - 1) + 5(1 - y 2 )] 2 

+ (a - b)[(a -b)(d-S)- m 2 b + Qp] 
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Table 1. Equations of state, used in Cook et al. 119941 
EOS Description 

A Reid soft core, adapted to nuclear matter 

AU Argonne V14 + Urbana VII 
FPS Urbana V14 + Three Nuclei Interaction 
L Nuclear attraction due to scalar exchange 

M Nuclear attraction due to pion exchange 
(Tensor interaction) 



y 2 f}{A(2Kmbx + 2m 2 b - Qp) 



(1 

[K 2 (x 2 -l)+S(l-y 2 )] + (l-y 2 ) 
{(a - b)(m 2 b 2 - 45d) - (4nmx + 2m 2 
-Q 2 )[(a - b)(d-8) - m 2 b + Qp]}). 



(17) 



Mathematically, the five parameters of this solution are 
independent. To make it applicable for astrophysical pur- 
poses and to describe astrophysical configurations, one needs 
to find constraints for the choice of and the relations be- 
tween some of these parameters. Those can only be found 
by matching the exterior solution with existing models for 
the interior of neutron stars. 



3 THE INTERIOR OF NEUTRON STARS 

Until now, the knowledge of the composition of the inte- 
rior of neutron stars is very poor. A large zoo of differ- 
ent EOS exists to describe all possible ingredients of the 
star and the interactions of these particles. Before the right 
EOS can be singled out by observations, i.e. determining 
the mass-radius-relation by e.g. measuring the cooling curve 
of neutron stars, the basic structure equations for rotating, 
axisymmetric stars - similar to the Oppenheimer-Volkoff- 
equations for spherical stars - have to be solved numerically 
for each model. Cook et al. 119941 performed these calcu- 
lations for 14 different EOS with a modified variation of 
the KEH-Code (Komats u et al. 198911 and tabulated several 
quantities for five of them representing the whole range of 
stiffness (A, AU, FPS, L, M; see table □ Cook et al. itTTOll 
and references therein). For further informations on differ- 
ent numerical methods see e.g. Ansorg et al. (2002} and 
references therein. For each EOS these authors computed 
models for three star sequences: first the normal sequence 
(NS) of stars whose masses are below the maximum mass 
of non-rotating stars, second the maximum mass normal 
sequence (MM) of stars with the maximum mass of non- 
rotating stars and third the supramassive sequence (SM) 
of stars whose masses exceed the maximum mass of non- 
rotating stars. These numerical values can be compared with 
the analytically derived ones to adjust the free parameters. 
The following fitting procedure was performed for all five 
equations of state and all three sequences. 

There exists a newer EOS including more effects of nu- 
clear and particle physics (Glendenning & Weber ( 200^ and 
references therein) , but unfortunately the analog analysis of 
Cook et al. (TjJSjJ has not yet been done. Perhaps this mod- 
ern EOS describes the interior of neutron stars in a better 
way than the older EOS considered by Cook et al., then the 
same fitting procedure performed with this EOS would give 



more physically meaningful results. More observations are 
needed to clear this point. 



4 FITTING THE FREE PARAMETERS 
4.1 The gravitational redshift 

One of the quantities tabulated in Cook et al. £1994) is the 
gravitational redshift photons experience being emitted for- 
ward (Zf) and backward (Zb) with respect to the rotation 
direction at the equator. The same can be easily calculated 
for this analytical solution. For photons which are emitted 
at the equator in forward (backward) direction with respect 
to the rotational direction the four-momentum has the form 

P * = const, x [ % + ^L_ € m )] (18) 



with the Killing- vectors and £f^v The frequency is cal- 
culated by 



where the fluid four-velocity on the equator is 

,-1/2 



(19) 



(20) 



and the fluid velocity measured by a ZAMO (zero angular 
momentum observer) 



v = ^/2wQ. -Q? (w 2 - p 2 /f 2 ). 

The frequency observed at infinity is 

Woo =Pni" t y 

Therefore the redshift which is looked for has the form 



(21) 
(22) 



(±2/pv / l-« 2 -2/ 3/2 (-l + v 2 )(TP + fw) 
+V 7 ! - v 2 y / 4f 2 p 2 v 2 - 4/ 4 (-l + v 2 )w 2 ^j I 



2 V 2 _ 4/4(-l + «2)w 2 

(2f 3 / 2 (-l + v 2 )( TP + fw)). 

A suitable combination of these two redshifts is 

1 



Z := (1 + Z f ) (1 + Z b ) 



f(Pe 



(23) 



(24) 



with the metric function / taken at the equator of the neu- 
tron star. To be able to compare the analytical function 



Zanal. = Z anaL (p, Z = 0; 771, a, 6, Q, fl) 



(25) 



with the numerical results one has to match the parame- 
ters. Cook et al. I|1994jl considered only neutron star models 
without electric charge and magnetic field, therefore one can 
set 



Q = 



and 



0. 



(26) 



Further it is useful to introduce the specific angular momen- 
tum j defined by 



J = 



./ 



a 
m 



m j. 



(27) 
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Because of the fact that Cook et al. (119941) made their calcu- 
lations in Boyer-Lindquist-like quasi-spherical coordinates r 
and 9 one has to transform the coordinates as 

r = | (V ' p 2 + (z + re) 2 + V P 2 + k) 2 ) + m 

i- (Vp 2 + (^ + ^) 2 - vV + (2 - *) 2 ) (28) 



u — arccos 

The connection of the distance coordinates in the equatorial 
plane (z — 0) is then 

m in 

(29) 



— (m + y/ pi + K 2 ^J 



= \J m 2 (r» — l) 2 



with the equatorial stellar radius Combining all these 
equations one gets 



b = b[Z(j),r(j);m(j),j] = b(j). 
It follows from the numerical data that 
M* = M„(j) = M,(-j) 

fl. = R*(j) = R*(-j) 



(30) 



(31) 



so b(j) also has to be symmetric with respect to the rota- 
tion direction. As we wanted to fit the parameter b in the 
whole range [-1,1] for j we had to include the Schwarzschild 
limiting case for j = 0. The solution of Manko et al. (2000 ) 
reduces to Kerr by setting 



(32) 



and to Schwarzschild with b — im. Therefore, we were es- 
pecially interested in imaginary values of b. This poses no 
problem as b is no physical observable. Moreover, it can be 
shown that the equation of the marginally-stable orbit (next 
section) has solutions for small j only with imaginary b. 



4.2 The marginally-stable orbit — an independent 
test 

The radius of the marginally-stable orbit is another quan- 
tity given in the tables of Cook et al. 1199411 . This provides 
the possibility to perform an independent test of the fits 
done above. The following calculation gives the analytic ex- 
pression for this radius. Consider a particle moving in the 
equatorial plane. With the three constants of motion E, L 
and K = 1 related to the energy, the angular momentum 
and the mass of the particle the equations of motion are 
fully integrable. The equation of the radial coordinate p can 
be interpreted as energy equation with an effective potential 



V(p)= e —f = ?--±(L-E W f-l. 
f TP 

For circular orbits the conditions 

v 0>) = y - ( L - E ») 2 - 1 

dV(p) _ 2/(L-Ew) 2 _ (L-Em) 2 /, p 
dp p 3 p' 2 



E 2 /, p | 2Ef(L-Ew)w,„ 



(33) 



= (34) 



(35) 



have to hold from which one can calculate 
E = v7 



L = E (w + X) 



(36) 

(37) 

p [-™,p ; 2 - yVp g + f, P p (2 7 =7Jp) ] rtsa 

X = / (2 /-/.„) ' (38) 

With the condition for a marginally-stable circular orbit 

d 2 V(p) _ -6f(L-Ew) 2 4(L-Ew) 2 /, P 
dp 2 p 4 p 3 

8E/(L-Ew)™, p 4E(L-Ew)/, p tt, p 



+ 



2E 2 / 2 2E 2 /^ 2 p E 2 /, p 

f 3 p 2 f 2 



(L-E</, w 2E/(L-Ew)w pp 

(39) 

and with 1361 . 137D and 1381 1 the equation for the radius of 
the marginally-stable orbit is 



(w, p w, pp f p (2/ - f, p p) + w 2 p f 4 [2f 2 + (-/» 

+f, PP f) P 2 ] + w, p f V™? P f 4 + f,pPW-f,pP) 
[2/ 2 + 2 f 2 p p 2 ~fp (4/, p + /, pp p)] + p(2f- f, p p) 
W,p /•' - 4/ 2 p fp + f% p 2 + f [f, pp p 

-™, PP f V™ 2 P f 4 + f,pp(2f-f. P p)]}) I 

[f P 2 {w 2 p f + 3/, p / p - f 2 p p 2 

-f [2 + w, p y/w 2 p f* + Tp P (2/ - f, P p)]}) = 0. 

(40) 

Again, one has to set z = Q — p — 0, then a function of 
m, r, j and 6 remains. To test the fitted parameters we used 
the m- and 6- values of the tables |3] and and made implicit 
plots in the (rj)-plane. 

5 THE FITS AND THE TESTS 

Fitting the numerically calculated masses in Cook et al. 
11994H with 



m(j) = m + mij 2 + m 2 j 4 + m 3 j 6 
and choosing the ansatz 



Ai a 2 ~m 2 + A 2 



+ A 3 



(41) 



(42) 



on finds the fitting parameters as listed in table |5] and ta- 
bic This ansatz is chosen to ensure 6 = im for j = 0, to 
achieve b(j) — b(—j) and to be sure that all terms have the 
right units. The fit parameters in 14H are in units of solar 
masses and the those in t|12t are dimensionless. Although b 
is imaginary the parameter re stays pure real. Therefore this 
solution is subextreme. Pure imaginary values of re represent 
hyperextreme spacetimes which could be interpreted as gen- 
erated by relativistic discs. The numerical data fall on the 
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Table 2. Dimcnsionless fitting parameters of the fit ansatz 



EOS A(NS) 



/ 

b = ./Ai a 2 - 
V 


m 2 + a 2 ^ 


+ A 3 ^| 




EOS (Seq.) 


Ai 


A 2 


A 3 


A(iNbj 


— (Ml /43 


—0.099 /8 


— 2.02 ( 92 


A /1\ /TA T\ 

A (MM ) 


—2.9820/ 


1 o noo rn 
— 12.02OD0 


14.2/300 


A /Q A 

A(bM j 


U.0lD04 


— 25. 00350 


30.50040 


AU(NS) 


-8.35758 


-7.59210 


11.33570 


AU(MM) 


-2.10796 


-7.94790 


8.71811 


AU(SM) 


0.83801 


-18.47420 


19.58550 


FPS(NS) 


-8.06897 


-20.34330 


40.12250 


FPS(MM) 


-3.94624 


-6.15936 


3.74551 


FPS(SM) 


-0.44467 


-21.45520 


25.25300 


L(NS) 


-19.75040 


6.10370 


-6.68697 


L(MM) 


-3.59342 


-7.91550 


7.25291 


L(SM) 


-0.10791 


-20.44890 


22.08320 


M(NS) 


-23.07410 


20.14150 


-37.30740 


M(MM) 


-6.50839 


-47.76620 


87.45900 


M(SM) 


11.14010 


-144.35100 


221.00000 


Kerr 


1 









Table 3. Fitting parameters in units of solar masses of the fit 
ansatz m(j) = mo + mij 2 + m,2j' 4 + m^j® 



EOS (Seq.) 


mo 


mi 


m,2 


m,3 


A(NS) 


1.39999 


0.08294 


-0.01655 


-0.00135 


A(MM) 


1.65510 


0.14606 


-0.02487 


-0.00866 


A(SM) 


1.73043 


0.20082 


-0.08822 


0.03529 


AU(NS) 


1.40004 


0.07455 


-0.01832 


0.00962 


AU(MM) 


2.13344 


0.21933 


0.03005 


-0.03461 


AU(SM) 


2.21414 


0.32712 


-0.12039 


0.07834 


FPS(NS) 


1.40003 


0.07255 


-0.02813 


0.02114 


FPS(MM) 


1.79945 


0.16119 


-0.04416 


0.01113 


FPS(SM) 


1.87636 


0.21989 


-0.11459 


0.05728 


L(NS) 


1.40000 


0.04002 


-0.00943 


0.00292 


L(MM) 


2.69993 


0.23158 


0.00001 


-0.03842 


L(SM) 


2.84355 


0.36050 


-0.19396 


0.11556 


M(NS) 


1.39999 


0.03695 


-0.01482 


0.00563 


M(MM) 


1.80458 


0.14357 


-0.27871 


0.31632 


M(SM) 


1.87629 


0.38213 


-0.91656 


0.89839 



theoretical curves with an accuracy of 0.6 - 1.8 % (see fig- 
ures ^^J. The small deviation is a consequence of rounding 
and numerical errors during the two fit processes of the mass 
and of b which could not be improved by taking higher orders 
into account. Because of the symmetry Z(j) = Z(—j) only 
the negative j'-values are plotted. Although this procedure 
was done for all sequences, we show only plots of the normal 
sequences, but the error range is valid for all sequences. 

The test of these parameters with the radius of the 
marginally-stable orbit shows deviations between the nu- 
merical data and the analytical function of the order 4 - 8 %, 
in the case of the normal sequence of the EOS L 15 % (Figs. 
EE)- In these fi gures the same function is plotted for the 
Kerr-metric (IMiller fc Lamb 1 996 ) for comparison. The ra- 
dius of the marginally-stable orbit of the normal sequence 
of the EOS M is always smaller than the stellar radius, so 
the corresponding plot is omitted. 
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Figure 1. gravitational redshift Z vs. rotation parameter j for 
the normal sequence of EOS A 
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Figure 2. gravitational redshift Z vs. rotation parameter j for 
the normal sequence of EOS AU 



6 SYMMETRY WITH RESPECT TO THE 
EQUATORIAL PLANE 

The ansatz of an imaginary parameter b was a direct con- 
sequence of the requirement that Schwarzschild is a limit- 
ing case for vanishing rotation. This choice causes a prob- 
lem. Following the formalism of deriving exact solutions 
of Einstein's equations the solution is symmetric with re- 
spect to the equatorial plane if certain multipole moments of 
the mass-, angular momentum-, electric and magnetic field- 
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Figure 3. gravitational redshift Z vs. rotation parameter j for 
the normal sequence of EOS FPS 
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Figure 4. gravitational redshift Z vs. rotation parameter j for 
the normal sequence of EOS L 
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Figure 7. radius of the marginally-stable orbit (in units of m) 
vs. rotation parameter j for the normal sequence of EOS AU 
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Figure 5. gravitational redshift Z vs. rotation parameter j for 
the normal sequence of EOS M 
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Figure 8. radius of the marginally-stable orbit (in units of m) 
vs. rotation parameter j for the normal sequence of EOS FPS 



distributions vanish (Manko ct al. 1994 1. This is obtained if 
the five parameters are real (IManko et al. 200 1 . Because of 
the imaginary parameter b, terms of the multipole moments 
change their role from mass to angular momentum multi- 
pole moments or from electric to magnetic field multipole 
moments and vice versa, and destroy this symmetry. One 
has to test if the physical gravitational and gravitomagnetic 
forces are still symmetric. Using the formalism of the 3+1- 
split (|Th orne et al. 1986 1 the four-dimensional space-time 
can be brought into the form 
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Figure 6. radius of the marginally-stable orbit (in units of m) 
vs. rotation parameter j for the normal sequence of EOS A 



ds 2 = a 2 dt 2 - h ik (dx l + $ dt) (dx k + f3 k di). (43) 

To perform the following calculations it is convenient to in- 
troduce a "star-fixed" coordinate system, at best in Boyer- 
Lindquist-like quasi-spherical coordinates. After comparison 
of 1431 1 with the metric Q in these coordinates 
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Figure 9. radius of the marginally-stable orbit (in units of m) 
vs. rotation parameter j for the normal sequence of EOS L 
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+[(r-m) 2 -k 2 ] sin 6 2 d<ff 



(44) 



- for the transformation of the coordinates p and z to r and 
6 see eq. 1281 - it follows 

2 _ / [(r - m) 2 - k 2 ] sin 6 2 



[(r — m) 2 — k?\ sin6> 2 — f 2 w 2 
f 2 w 



[(r — my — sm f 
/3 r = f3 e = 
e 27 [fr - 



f 2 w 2 



/ [(r — m) 2 — k 2 ] 
e 27 [(r-m) 2 -k 2 cos<9 2 ] 
7 

[(r - m) 2 - k 2 ] sme 2 -/ 2 ^ 2 



/ 







(45) 



The gravitational and gravitomagnetic force can be calcu- 
lated with 

1 u 1 a 



9 



(46) 



which means 



H. 



64> 



H 



tie 



H 




with the orthonormal basis vectors 
1 d 



(47) 




(48) 



Inserting all functions in 1471 . one finds that the transfor- 
mation r,9 = ^ — ■# — > r, 6 = 5+1? does not change the coef- 
ficient of e.f and the components H.^ and H^-, but changes 
the sign of the coefficient of eg and the components Hg^ 
and Hig. This shows that the physical forces are symmetric 
with respect to the equatorial plane. 



7 THE ERGOREGION 

The existence of the ergoregion is a property of stationary, 
axisymmetric space-times which is defined as the region be- 
tween an event horizon (if existing) and the stationary limit 
surface. In stationary, axisymmetric space-times the compo- 
nents (? M „ of the metric tensor are independent of the cyclic 
variables t and 4>. Therefore two Killing vectors 



^ = i und *<« = w 

are existing. With the angular velocity 

d(j} d(f> I dr 



(49) 



u 
a* 



dt dt/dr 

the four-velocity of a stationary observer takes the form 

u" - uH d i n 9 )- *w +ng w 
- U( st + V ^, + ^,1 



(50) 



(0) 



with which the condition for SI 



(51) 



W 2_5!i w = _^i (52 ) 



si_ < si < si+ si± = u ± 



is introduced. There is a limit where no static observer with 
S7 = can exist. The surface of this stationary limit can be 
calculated from 



gu = o 



/ = o. 



(53) 



The implicit plots of this surface reveal new features 
which are unknown for the Kerr-solution. In contrast to the 
Kerr-metric, the ergoregion inflates with increasing rotation 
parameter j for neutron stars in their normal sequence (their 
mass is below the maximum of a non-rotating star), while 
the EOS and the rest mass of the star are kept fixed (Fig. 

rrm . 

Additionally the ergoregion consists beside the main 
body of two "bubbles" lying at the poles of the ergoregion 
whose existence is completely unknown in the Kerr-metric. 
As the rotation speed increases, the bottle neck connecting 
the bubbles with the main ergoregion attenuates until, at a 
critical , these bubbles peel away. 

Another possible effect of rapid rotation is a constriction 
of the ergoregion in the equatorial plane which also increases 
for higher rotation. At another critical rotation value jc°S str 
two separated regions appear (Fig. lllH . This takes place only 
on neutron stars with stiffer EOS FPS, L and M in the 
normal sequence (notation see Cook et al. |1994|). 

Until now, the magnetic dipole moment has been ne- 
glected. The surprising result of an increased non- vanishing 
magnetic field is the fact that the main ergoregion is almost 
not affected fFig. 11211 . while the bubbles grow in size and the 
constriction moves to larger values of z coordinate (Tig. liaH . 
Because of the fact that observed neutron stars have mag- 
netic fields in the range of 10 12 - 10 15 gauss, the parameter 
p should be smaller than 0.1. Therefore, the effect of a non- 
vanishing magnetic field is relatively small. Neutron stars 
are flattened due to their rotation and oblate; their equato- 
rial radius is larger than the distance between their poles, 
whereas the ergoregion is prolate. Three different cases are 
then possible: 
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Figure 10. Inflation of the ergoregion with increasing rotation 
parameter j 
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Figure 11. Constriction of the ergoregion in the equatorial plane 
and splitting into two main ergoregions at a critical rotation 




0) : 
star 

(ii) 



> r* , the ergoregion is completely outside the 



.erg 
eq 



< rZ 



but 



c rq 
pole 



> r 



pole i 



the ergoregion exceeds 



the stellar surface only in the polar region 



(iii) 



< r* eq and r po f e 



< r 



pletely inside the star. 



pole ' 



the ergoregion is com- 



Table|l]shows which case occurs for which EOS if j is smaller 
than the j at mass-shedding for uniform rotation. For all 
EOS in the normal sequence the ergoregion is inside the 
region for slow rotation, where the stellar surface lies in the 
models of the interior. Increasing j, the ergoregion moves 
outside in the polar region after passing a critical j^f t . Only 
stars with the EOS AU can be accelerated to this critical 
value, all other stars become unstable due to mass-shedding 
at smaller j. The critical j decreases with increasing stellar 
mass and increases with increasing stiffness. 
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Figure 12. Invariance of the main ergoregion to an increasing 
magnetic dipole moment fi 
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Figure 13. Effects of an increasing magnetic dipole moment p. 
to the polar bubbles 



Table 4. Position of the ergoregion in comparison to the stellar 
surface (Q = 0, p = 0) 
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A normal sequence 


0.3 


(iii) 




0.6 


(hi) 


A maximum mass normal sequence 


0.3 
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0.6 


(ii) 


A supramassive sequence 


0.6 


(") 


AU normal sequence 


0.3 


(iii) 




0.6 


(ii) 


AU maximum mass normal sequence 


0.3 


(ii) 




0.6 


(ii) 
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0.6 
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FPS normal sequence 
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0.6 


(iii) 


FPS maximum mass normal sequence 


0.3 


(ii) 
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(iii) 


M maximum mass normal sequence 
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(iii) 


M supramassive sequence 
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(ii) 



8 SINGULARITIES OF THE GLOBAL 
SOLUTION 

Singularities are given if the curvature tensor diverges. This 
divergence occurs if the denominator of the Ernst-potentials 
IKrnst 196811 has to vanish JManko et al. 199911 : 

A + 2mB = (54) 

with 

A = 4[(kV -6y 2 ) 2 -d 2 -iK S xy{a-b)(x 2 - 1)] 
-(1 - y 2 )[{a -b)(d-8)- m 2 b + Qp] 
[(a - b)(y 2 + 1) + iinxy] 
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Figure 14. Splitting of the point singularity (j = 0) into two ring 
singularities (j ^ 0); only the singularity below the equatorial 
plane is shown (z < 0) due to clarity (dashed line) 
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Figure 15. Same as Fig. El but also merging of the two ring 
singularities to one single ring singularity (p ^ 0, z = 0) at \j\ m 
0.9 

B = kx{2k 2 (x 2 - 1) + [b(a - b) + 2<5](1 - y 2 )} 
+iy{2K 2 b(x 2 — 1) — [k 2 (a — b) — m 2 b + Q\x 
-2a5](l-y 2 )}. (55) 

Equation 1541 in fact consists of two equations, because A 
and B are complex and both, the real and imaginary part, 
have to vanish simultaneously. These two equations give the 
coordinates x-y, or p-z respectively, of the singularities. 

These were computed numerically for all EOS and se- 
quences as a function of the rotation parameter j. If j = 
only one point singularity arises in the origin as expected 
for the Schwarzschild limiting case. With increasing rota- 
tion this singularity splits into two ring singularities whose 
radius increases and which move away from the equatorial 
plane (Fig. 1141 . In a few cases the singularities merge again 
to one single ring singularity at a critical (Fig. ll5ll . but 
the data of Cook et al. Q994 ) suggest that stars become un- 
stable before reaching this rapid rotation. Again the results 
shown apply only to non-magnetized stars. In the presence of 
magnetic fields the two ring singularities are present (p ^ 0, 
2 7^ 0) even for vanishing rotation. The Schwarzschild space- 
time is of course no limiting case anymore. Increasing j, the 
singularities approach each other, merge at a critical rota- 
tion j™?t UB - but not in the origin - (sign change of p while 
and seem to change their roles (sign change of z) at 
jcrtt- The singularity lying above the equatorial plane for 
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Figure 16. Radius of the ring singularities; two ring singularities 
for vanishing rotation (p ^ 0), merging at a critical rotation (sign 
change of p) 
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Figure 18. Ring currents as sources of the exterior fields 
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Figure 17. Distance of the ring singularities to the equatorial 
plane; moving of the singularities from above to below the equa- 
torial plane (sign change of z) 

slow rotation moves after the merging below the equatorial 
plane (Fig. fEl EJ ■ 



9 HORIZONS 

The horizon is the null surface spanned by the tangential 
vectors dt and d$. The resulting condition is 

gttg^-glt = -p 2 = 0- (56) 

p plays in this coordinate system the role of the Weyl radius. 

The singularities are located at p > 0. Therefore the 
solution hurts the "cosmic censorship hypothesis". But be- 
cause of the fact that all singularities are within the stellar 
surface this is irrelevant. 



10 CONCLUSIONS 

The same procedure as shown in section 0] was already per- 
formed - in the opposite order - by Sibgatullin & Sunyaev 
( 1998), but only with the normal sequences of EOS A, EOS 
AU and EOS FPS and with an older analytic solution found 
by Manko et al. 119941 which did not permit a rational func- 
tion representation and could not be written in a concise 
form. These three EOS are of lower stiffness. These authors 



arrived at the conclusion that it is sufficient to include only 
the mass quadrupole moment and no higher multipole mo- 
ments - a result we could not confirm for all, especially not 
for the stiffest EOS L and M. The fact that the accuracy is 
high in the case of the redshift, but much poorer in describ- 
ing the radius of the marginally-stable-orbit which is more 
sensitive to higher multipole moments of the space-time (see 
comparison to Kerr in the figures) , seems to show that higher 
multipole moments have to be included, especially for rapid 
rotation and stiff er EOS, to improve the form of the ex- 
terior gravitational field of neutron stars. But until a new 
exterior solution is found, this new self-consistent relativis- 
tic model of rapidly rotating neutron stars, we present here, 
is a great improvement in comparison to the Schwarzschild- 
and Kerr-solutions sometimes used to describe also rapidly 
rotating stars. 

It is clear that a vacuum solution is not able to de- 
scribe the interior of the star. For each EOS the parame- 
ters have to be adjusted to match both regions. Therefore 
this solution cannot claim to be global and the presence 
of (naked) singularities which are always inside the stellar 
surface is irrelevant for astrophysics. However they help to 
interprete the exterior field as a space-time created by a 
configuration with two ring currents (Fig. 1180 . These mass- 
and angular momentum currents arise as singularities in the 
curvature tensor. Coupled by the Einstein equations to the 
energy-momentum-tensor of the electromagnetic fields these 
sources should also appear in the electric and magnetic field- 
distributions. Because of charge neutrality the ring currents 
must have different rotation directions and therefore create 
also weak electromagnetic quadrupole fields. As discussed in 
section [7] it seems to be common that the ergoregion is not 
completely enclosed by the stellar surface. The fact that it 
leak from the interior of the star only in the polar regions 
could have very interesting consequences for the acceleration 
of particles and propagation of light. 
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